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Resumen 


Cortés-Martínez, F., Treviño-Cansino, A., Alcorta-García, M. 
A., Sáenz-López, A., £ González-Barrios, J. L. (marzo-abril, 
2015). Optimización en el diseño de lagunas de estabilización 
con programación no lineal. Tecnología y Ciencias del Agua, 
6(2), 85-100. 


El estudio presenta un modelo matemático de optimización 
para el diseño de sistemas lagunares integrado por dos 
lagunas: facultativa y de maduración. El modelo tiene 
como función objetivo el costo y considera cuatro variables 
de decisión: tiempo de retención hidráulico y número 
de mamparas en la laguna facultativa; en la laguna de 
maduración, también número de mamparas y tiempo de 
retención. Se consideran como restricciones los coliformes 
fecales y la materia orgánica. Dicho modelo resultó ser 
no lineal, ya que las relaciones entre las variables no son 
proporcionales. Se diseñó un sistema lagunar utilizando la 
metodología tradicional adoptada por la Comisión Nacional 
del Agua para México, luego se aplicó el modelo matemático. 
Se lleva a cabo una comparación de los resultados, los cuales 
indican una reducción del tiempo de retención de 14.16% y 
una disminución del costo de 12.04%. El resultado anterior es 
importante, ya que la principal desventaja de estos sistemas 
es el requerimiento de terreno. Se incluye, además, el análisis 
de sensibilidad de la función objetivo y las restricciones 
consideradas. El citado análisis es sensible a la variación 
de los parámetros. Ambos estudios cumplen perfectamente 
con las condiciones de calidad del agua tratada que indica 
la normatividad para el vertido a los cuerpos receptores. Se 
recomienda llevar a cabo otros estudios de optimización, 
considerando diferentes configuraciones. 


Palabras clave: sistemas lagunares, coliformes fecales, 
tiempo de retención, materia orgánica, costo de construcción, 
modelo matemático de optimización. 


Abstract 


Cortés-Martínez, F., Treviño-Cansino, A., Alcorta-García, M. A., 
Sáenz-López, A., € González-Barrios, J. L. (March-April, 2015). 
Design Optimization for Stabilization Ponds Using Non-Linear 
Programming. Water Technology and Sciences (in Spanish), 
6(2), 85-100. 


The article herein presents a mathematical optimization model 
to design lagoon systems containing two ponds: facultative and 
maturation. The objective function of the model was the cost. Four 
decision variables were considered: hydraulic retention time and 
number of screens in the facultative pond; in the maturation pond, 
also number ofscreens and retention time. Fecal coliform and organic 
matter were considered as restrictions. The resulting model was non- 
linear since the relations among the variables were not proportional. 
A pond system was designed using the traditional methodology 
adopted by the National Water Commission for Mexico. The 
mathematical model was then applied. A comparison of the results 
indicated a reduction in retention time of 14.16% and a decrease in 
cost of 12.04%. This finding is important given that the primary 
disadvantage of these systems is the need for land. In addition, a 
sensitivity analysis of the objective function and the restrictions is 
included. This analysis is sensitive to changes in the parameters. 
Both studies fully meet the quality standards stipulated by norms 
for treated water discharged into receptor bodies. Additional 
optimization studies that include different configurations are 
recommended. 


Keywords: Pond systems, fecal coliform, retention time, organic 
matter, construction cost, mathematical model optimization. 
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Introducción 


Las lagunas de estabilización tienen tres objeti- 
vos principales: a) la eliminación de coliformes 
fecales; b) la remoción de la materia orgánica, 
también llamada demanda bioquímica de 
oxígeno (DBO), y c) de nutrientes: nitrógeno y 
fósforo (CNA e IMTA, 2007a). De acuerdo con 
Senzia, Mayo, Mbwette, Katima y Jorgensen 
(2002), Agunwamba, Ochonogar y Ukpong 
(2003), Mara (2004), Abbas, Nasr y Seif (2006), 
Hamzeh y Ponce (2007), y Naddatfi et al. (2009), 
los estanques de estabilización se sugieren en 
países con climas tropicales, pues las condicio- 
nes ambientales incrementan la eficiencia en la 
remoción de contaminantes. 

Los sistemas lagunares se pueden clasificar 
en relación con la presencia de oxígeno: anae- 
robias, facultativas y de maduración (Metcalf éz 
Eddy, Inc., 1991): 


e  Anaerobias. Las bacterias presentes no 
requieren oxígeno disuelto necesario para 
la descomposición de la materia orgánica; 
al proceso descrito también se le conoce 
como digestión anaerobia (Rolim, 2000). 

e  Facultativas. El mecanismo se lleva a cabo 
en el estrato superior, es decir, el comen- 
salismo de bacterias aerobias y algas. La 
materia orgánica es descompuesta por 
las bacterias heterotróficas, lo anterior ge- 
nera compuestos inorgánicos. El oxígeno 
necesario para realizar la simbiosis es 
suministrado principalmente por el proce- 
so de fotosíntesis (Mara, 2004). 

e De maduración. Reciben el efluente de las 
lagunas facultativas y se utilizan para pulir 
el efluente de patógenos de acuerdo con 
la calidad requerida. Sólo existe una zona: 
aerobia (CNA e IMTA, 20076). 


Deflectores 
Shilton y Mara (2005), Abbas et al. (2006), Cor- 


tés, Treviño, Luévanos, Luévanos y Uranga 
(2014a), y Cortés et al. (2014b) recomiendan el 


uso de deflectores en los sistemas lagunares, 
puesto que se mejoran en forma importante las 
condiciones hidráulicas y se favorece la reduc- 
ción de las áreas muertas. 

A la fecha se han realizado estudios en labo- 
ratorio para aumentar la eficiencia hidráulica 
y la eliminación de contaminantes, conside- 
rando deflectores. Killani y Ogunrombi (1984); 
Pedahzur, Nasser, Dor, Fattal y Shuval (1993); 
Muttamara y Puetpaiboon (1997); Zanotelli, 
Medri, Belli-Filho, Perdomo y Costa (2002); 
Shilton y Harrison (2003a); Sperling, Cherni- 
charo, Soares y Zerbini (2003); Shilton y Mara 
(2005), y Abbas et al. (2006) publicaron que se 
obtiene mayor eficiencia hidráulica y mejor 
tratamiento del agua residual con mamparas 
a 70% de la longitud. Luego Banda (2007), y 
Winfrey, Stronsnider, Nairn y Strevett (2010) 
concluyeron que se obtienen mejores resulta- 
dos al utilizar mayor número de deflectores. 


De los modelos matemáticos 


Killani y Ogunrombi (1984) recomendaron 
el uso de la investigación de operaciones en 
lagunas de estabilización, con el propósito de 
obtener una optimización de recursos econó- 
micos, es decir, determinar el costo óptimo del 
sistema. Nelder y Mean (1965); Fonseca y Fle- 
ming (1993, 1995); Oke y Otun (2001); Bracho, 
Lloyd y Aldana (2006); Winfrey et al. (2010), y 
Olukanni y Ducoste (2011) aplicaron la progra- 
mación lineal en el diseño de sistemas lagu- 
nares; concluyeron que sí era posible mejorar 
el diseño mediante la optimización, a manera 
de incrementar la eficiencia del tratamiento y 
minimizar costos. Finalmente, Sah, Rousseau y 
Hooijmans (2012) realizaron estudios compara- 
tivos de los modelos existentes y concluyeron 
que era necesario proponer un modelo integral 
de optimización, donde se consideren todas las 
variables que intervienen en el tratamiento del 
agua con estanques de estabilización. 

De acuerdo con la revisión bibliográfica, en 
México existen pocos modelos matemáticos de 
optimización aplicados a los estanques de es- 
tabilización; debido a lo anterior, se plantearon 
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los siguientes objetivos: a) realizar un análisis 
matemático para determinar el modelo de 
optimización con dos lagunas; b) diseñar un 
sistema lagunar con la metodología tradicio- 
nal, incluyendo el costo; c) aplicar el modelo 
matemático de optimización utilizando el sis- 
tema Solver de Excel, y d) comparar los resulta- 
dos entre los dos análisis y definir las ventajas 
económicas de ambos estudios. 


Programación no lineal 


El Gradiente Reducido Generalizado (GRG) es 
el algoritmo utilizado por el sistema Solver de 
Excel. El GRG inicia con una solución conocida, 
es decir, calculada anteriormente. El objetivo es 
encontrar una nueva propuesta, pero mejora- 
da. El algoritmo revisa la región factible de la 
solución conocida y determina una nueva solu- 
ción, hasta que ya no es posible optimizar más 
el resultado. La herramienta actúa en un grupo 
de celdas llamadas variables de decisión, que 
se incluyen en la fórmula de la función objetivo 
y restricciones. Algunas de las características 
de los problemas de programación en donde 
es posible utilizar el método de GRG son las 
siguientes: la función objetivo puede ser lineal 
o no lineal, admite restricciones de igualdad o 
desigualdad y el punto de inicio para aplicar 
el GRG puede ser factible o no (Muramatsu, 
2011). 


Análisis de sensibilidad 


Según Anderson, Sweeney y Williams (2004), 
el diagrama de tornado es recomendable para 
determinar la sensibilidad del proyecto. El 
análisis consiste en modificar los rangos de va- 
lores considerados en el estudio, condición que 
proporciona información acerca de cómo afecta 
a la solución mejorada. El diagrama de tornado 
muestra la información por medio de barras, es 
decir, la barra más ancha indica el parámetro 
más sensible al cambio. 

El presente estudio es continuación del 
modelo matemático de optimización publicado 
por Olukanni y Ducoste (2011). La aportación 


del presente artículo es determinar un modelo 
matemático, incluyendo dos lagunas: facultati- 
va y de maduración. Se consideran como res- 
tricciones la demanda bioquímica de oxígeno y 
el número de coliformes fecales en el efluente. 

Consideraciones del modelo matemático 
publicado por Olukanni y Ducoste (2011): 
utilizan un modelo multiobjetivo; es decir, se 
maximiza la eliminación de contaminantes y se 
minimiza el costo; incluye además como res- 
tricción la concentración de coliformes fecales 
en el efluente, el número de mamparas y el área. 
Mientras que en el presente estudio se incluyen 
como limitaciones tanto los coliformes fecales 
como la materia orgánica, número de mam- 
paras y tiempo de retención. Otra diferencia 
importante es que el presente estudio incluye 
un solo modelo para dos lagunas, mientras que 
el de Olukamni considera una sola laguna. 

La novedad del presente trabajo respecto a 
otras publicaciones de los autores referencia- 
dos (Cortés, Treviño, Luévanos € Luévanos, 
2013; Cortés et al., 2014a; Cortés et al., 2014b), 
como ya se indicó, es que se considera un mo- 
delo matemático integrado por dos lagunas, 
con cuatro variables de decisión; mientras que 
en las publicaciones anteriores, el modelo fue 
definido para una sola laguna, con dos varia- 
bles controlables. Es importante indicar que el 
presente modelo matemático es único, es decir, 
de acuerdo con la revisión bibliográfica, no ha 
sido publicado anteriormente. 

El documento está organizado de la si- 
guiente forma: en la primera parte se describe 
la nomenclatura para el diseño de una laguna 
facultativa y de maduración, así como la me- 
todología tradicional de diseño. La segunda 
parte incluye el análisis realizado para definir 
la función objetivo preliminar, relaciones ma- 
temáticas, variables de decisión y restricciones 
consideradas. En la tercera parte se desarrolla 
un ejemplo de diseño con la metodología tra- 
dicional y luego la aplicación del modelo ma- 
temático propuesto; además se considera un 
cuadro comparativo con resultados y análisis 
de sensibilidad del modelo. En la cuarta parte 
se incluyen las conclusiones. 
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Materiales y métodos 


Para el diseño del sistema lagunar se consideró 
la metodología establecida para México por la 
Comisión Nacional del Agua (CNA) y el Insti- 
tuto Mexicano de Tecnología del Agua (IMTA). 
El manual consultado fue Paquetes tecnológicos 
para el tratamiento de excretas y aguas residuales 
en comunidades rurales de la Comisión Nacional del 
Agua (CNA e IMTA, 2007a). 


Nomenclatura para el diseño de lagunas 
de estabilización 


CO — = carga orgánica. 

Q, = caudal en el influente de la laguna 
(m? / día). 

DBO, = concentración de la demanda bio- 
química de oxígeno en el influente 
de la laguna en mg/l. 

1000  = factor de conversión. 

Av = carga orgánica superficial (kg/ha- 
día). 

E = temperatura mínima media men- 
sual del aire (*C). 

L, = concentración de DBO, en el in- 
fluente de la laguna (mg/1). 

Af = área de la laguna facultativa en (m?). 

Omed = caudal en el influente en (m? / día). 

£ = profundidad (m). 

V = volumen (m?). 

EL = tiempo medio de retención hidráuli- 


co, laguna facultativa (días). 


Ol = tiempo medio de retención hidráuli- 
co, laguna de maduración (días). 

XxX = relación entre longitud y anchura. 

Boro  = ancho promedio en (m). 

Loi.  = longitud promedio en (m). 

Bono = ancho superior en (m). 

L¿,,  = largo superior en (m). 

A = área superficial en (m?). 

Qe = caudal en el efluente de la laguna 
(m? / día). 

e = evaporación (mm / día). 

d = factor de dispersión adimensional. 

K, = coeficiente de reducción bacteriana 
(a). 


a = constante adimensional. 

Ne = coliformes fecales corregidos por 
evaporación en el efluente de la 
laguna (NMP/100 ml). 

Ni = coliformes fecales en el influente de 
la laguna (NMP/100 ml). 

Nf/No = número de coliformes fecales en el 
efluente (NMP /100 ml). 

Kf = constante de decaimiento de la 
DBO, (día”). 


DBOef = concentración de la DBO, en el 

efluente de la laguna (mg/1). 

DBOe = concentración de la DBO, en el 
efluente de la laguna corregida por 
evaporación (mg/!). 

= número de mamparas en la laguna 


MampF 
facultativa. 
Mampm Número de mamparas en la laguna 
de maduración. 
An = área de la laguna de maduración 


(m?). 


Diseño de laguna facultativa (flujo 
disperso. Método de Yánez) 


a) Carga orgánica: 
CO. =Q¡(DBO,)/1 000 (1) 
b) Carga superficial de diseño: 
As = 250(1.085)'=" Q) 
c) Área de la laguna facultativa: 
Af =10L,Qmed/hs (3) 
d) Volumen de la laguna: 
V =(AfAZ) (4) 
e) Tiempo medio de retención hidráulico: 
O,¿=V/Q, (5) 


f) Dimensionamiento. Relación largo ancho 
X=3: 
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Brant ¡E 6) 
E 0) 
g) Para el ancho y largo: 
Bóup = Brom + (2)Talud) (8) 
os = Lo, +(Z)(Talud) (9) 
h) Cálculo del área superficial: 
Aso = (Ba, XL.) (10) 
1) Gasto en el efluente: 
Qe=Q, -0.0014,,e (11) 


j) Remoción de coliformes fecales: consideran- 
do mamparas con longitud de 0.70 del largo de 
la laguna se tiene que: 


X = (L,,...) (0.70) (núm. de mamparas 
+ 1)/(B,,.,,)/ (núm. de mamparas +1) (12) 
ds e 7 (13) 
-0.26118 +0.25392(X) +1.0136(X) 
k) Coeficiente de reducción bacteriana: 
K, =0.841(1.075)2 (14) 


MAN 


1) Constante “a”. Se determina a través de la 


fórmula: 


a=,/1+4(K,0,d) (15) 
m) Coliformes fecales en el efluente de la lagu- 
na facultativa: 


Nf _ 4aexp 02 


j 16 
No (1+a) dl do 


n) Coliformes fecales corregidos por evapora- 
ción: 


Ne =(NeX(Q;) / Qe (17) 
0) Concentración de la DBO en el efluente de 
la laguna: 


Kf a Kf. 35 


(1.085) 5" Sd 


p) Concentración de la demanda bioquímica de 
oxígeno en el efluente de la laguna: 


Dios. 22 (19) 
KfO, +1 
q) Eficiencia de remoción de la DBO: 
ae (20) 
DBO, 
r) DBO corregida por evaporación: 
DBOe=(DBO'(Q,) / Qe (21) 


Laguna de maduración (método flujo 
disperso) 


1. Tiempo de retención hidráulico (O). Se 
analiza por tanteos, es decir, se propone un 
tiempo de retención hidráulico y se revisa 
la concentración de los coliformes fecales 
en el efluente; éstos deben ser iguales o 
menores de 1000 NMP/100 ml. 


2. Volumen de la laguna: 


V =(QNOy) (22) 
3. Área de la laguna: 
Vv 
Au=z (23) 


Para el resto del diseño se aplicaron las 
ecuaciones a partir del inciso f) hasta r), que se 
indican en la metodología de la laguna facul- 
tativa. 
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Función objetivo preliminar 


Los costos considerados para determinar la 
función objetivo son los siguientes: del terreno, 
$750.00 metro cuadrado; costo índice del bordo, 
$1 200.00 metro lineal; de la mampara, $500.00 
metro lineal. 

Según Cortés et al. (2014a; 2014b), la longi- 
tud de la mampara se considerará al 70% de 
la longitud: (0.7) Na Para definir el área 
de la laguna facultativa se observa el ancho 
superior más la longitud de cada uno de los 
lados, como se muestra en la figura 1. Para la 
longitud se tiene Bs + 12) E 12); de igual 
forma para el perímetro: 2 Ba, 3) +2 E 3). 
El mismo criterio se sigue para la laguna de 
maduración, considerando la información que 
se indica en la figura 2. Con los datos anteriores 
se determinó la expresión (24): 


Costo total =750|(B,,,,,-+12)(L,p +12) 


+(Losprr+ 10) (Lagar +3)] 
Hl 200|2(Ba,p,+ 3) +2(Lo,p,+3)+ (Boypar+ 3) 


+2 (1, +3)] 


SupM 


MampF""SupF MampM=SupM 


A A A A (24) 


Donde Bono =Brromy+3. Según la CNA y el 
IMTA (2007a; 2007b), esta expresión se refiere 
al ancho superior de la laguna facultativa más 
tres metros de la relación 2:1 considerada en los 
bordos (figura 1). La ecuación siguiente deter- 
mina la relación largo ancho: 


Bóup = ancho superior 


Profundidad 1.50 m 


Brrom = ancho promedio 


Figura 1. Corte transversal de la laguna facultativa. 


Bsup = ancho superior 


Profundidad 1.00 m 


Brrom = ancho promedio 


Figura 2. Corte transversal de la laguna de maduración. 
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Para calcular el ancho promedio de la 
laguna facultativa se realizan las siguientes 
operaciones: a) despejar el volumen de la 
ecuación (5) de la metodología tradicional; 
b) luego sustituir el volumen en la expresión 
(4), y c) finalmente la expresión se determina 
para el ancho del estanque, sustituyendo en la 
ecuación (6): 


Para calcular el ancho superior de la laguna 
se considera la relación de talud de los bordos; 
de esta forma se define la expresión (25): 


=_./0r*Qr 
Boupr SS 4.5 + 3 


(25) 
De forma similar, para el largo superior de 


la laguna facultativa se obtiene la expresión 
(26) (Cortés et al., 2014a, 2014b): 


22 43=,/20,*Q,+3 


Para la laguna de maduración con 1.0 metro 
de profundidad se define: 


Esupr 7 3 (26) 


B B 1 


SupM e SupF dl 


Sustituyendo la expresión (25) en 


B. .—1se obtiene: 


SupF 
= el 0,*0j 
Boro 0 45 + 2 


Para calcular el largo de la laguna de madu- 


B 


SupM 


(27) 


ración se considera la ecuación (7) de la meto- 
dología; luego, siguiendo el mismo criterio del 
ancho superficial de la laguna de maduración, 
se tiene: 


Prom + 2 
PromM 


Para el área se sustituye la ecuación (22) en 
la expresión (23) de la metodología; como ya 
se dijo, se considera 1.0 metro de profundidad 


para la laguna de maduración. El ancho pro- 
medio se obtiene restando la profundidad por 
la relación de inclinación del talud al ancho 
superficial. Finalmente se obtiene la expresión 
para calcular a longitud superior: 


(28) 


Sustituyendo las ecuaciones (25), (26), (27) y 
(28) en la (24), se determina la función objetivo, 
ecuación (29): 


Costo total =750|[/S:% +15)[/20,+Q, +15) 
22 +12) 
+1200|2//%2 +6]+2(/20,=0, +6) 
EE +5)+2 e) 
+500(0.7)|N nos (20,+0, +3) 


+ 


4.5 


| (29) 


Como se observa en la expresión (29), el mo- 
delo matemático resultó no lineal debido a las 
relaciones matemáticas que intervinieron para 
determinar la función objetivo. 

Las variables de decisión consideradas son 
Os, y ¿Nuramprs N mampm! 

El tiempo de retención hidráulico en ambas 
lagunas influye directamente en las dimensio- 
nes de éstas. El número de mamparas favorece 
la eficiencia del tratamiento de las aguas 
residuales (Winfrey et al., 2010). 

Restricciones: 


1<0,<30 


d>0 
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1<0,,s10 
0O<DBOe<75 


Ne =1000 


N 


MampF Y 


N 


Mampm = enteros 


Nuampr/ Op 30 

Se diseñó una planta de tratamiento de 
lagunas de estabilización para un poblado 
ejidal. La comunidad en cuestión se ubica en 
el municipio de Gómez Palacio, Durango, Mé- 
xico. Datos: número de habitantes al horizonte 
de proyecto 1 500; aportación = 154 1/hab /día 
(CNA e IMTA, 2007a); gasto de diseño: 231 m?/ 
día; temperatura promedio del mes más frío: 
11.8 “C (CNA e IMTA, 2007b); concentración de 
la materia orgánica (DBO,) en el influente =220 
mg/l; coliformes fecales 1.0 x 107 NMP/100 ml. 
Estos valores son considerados para aguas resi- 
duales domésticas (Metcalf €: Eddy, Inc., 1991). 


La evaporación es de 5 mm/ día. Considerar la 
relación largo ancho de tres para la laguna fa- 
cultativa; para la de maduración se considerará 
el mismo ancho que el estanque facultativo. 
Para determinar la calidad del agua residual 
tratada se incluyó como base lo que indica la 
Norma Oficial Mexicana NOM-001-ECOL-96 
(DOF, 1996). 


Resultados y discusión 


El cuadro 1 muestra los resultados de diseño 
del sistema de lagunas con el método tradicio- 
nal; la figura 3 indica las dimensiones deter- 
minadas; se observa que con 2 y 4 mamparas 
se cumple con las condiciones de calidad que 
indica la norma. 


Aplicación del modelo matemático para 
la optimización 


El cuadro 2 indica que el sistema encontró una 
solución en donde todas las restricciones y con- 
diciones de optimización fueron satisfechas. En 
el primer apartado del cuadro 2 se observa el 
comparativo de costos. El resultado fue 12.04% 
menor; en el mismo cuadro, en el apartado de 
variables de decisión, se indica que el sistema 


Cuadro 1. Resultados de diseño del sistema lagunar, método tradicional. 


Datos Laguna facultativa Laguna de maduración 
Qi Ni DBOi T O, NN nampr Os, Nod 
28il 10 000 000 220 11.8 AN 2 1.88 4 
Resultados de la laguna facultativa 
Xx d kb a Wa, E Qe Ne DBOe Área 
18.9 0.0516 0.4648 1.8628 S3/ 112.11 208.93 2 ES 43 4 414.03 
Resultados de la laguna de maduración 
Xx d kb a Wa, Lo. Qe Ne DBOe Área 
58.94 0.0167 0.4648 1.0287 38.37 12.80 206.47 998.09 32 491.13 
Área total | Costo, laguna facultativa | Costo, laguna de maduración | Costo total 
4 905.16 $5 238 328.77 $678 661.09 $5 916 989.86 
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Dirección de flujo 


Efluente DBO; = 32 mg/l 
Coliformes fecales = 998 NMP /100 ml 


Maduración 


Longitud de la mampara, 78.47 


39.37 m Laguna facultativa 


38.37 m 


112.1 m 


Influente DBO; = 220 mg/l 
Coliformes fecales = 1.0 x 107 NMP/100 ml 


12.8m 


Influente DBO; = 43 mg/l 
Coliformes fecales = 2 335 NMP/100 ml 


Figura 3. Dimensiones del sistema lagunar con el método tradicional. Arreglos de la laguna facultativa con dos mamparas 


y de la laguna de maduración con cuatro. 


redujo el tiempo de retención hidráulico de 
25.77 a 22.12 días: la diferencia es de 3.65 días, 
que representan 14.16%. Según la CNA e IMTA 
(2007a; 2007b), el tiempo de retención hidráuli- 
co influye de forma directa sobre el tamaño del 
estanque y, como consecuencia, en la necesidad 
del terreno. Los resultados del presente docu- 
mento coinciden con las aseveraciones de la 
CNA e IMTA (2007a; 2007b). Otro cambio impor- 
tante es que se determinaron cuatro mamparas 
en la laguna facultativa, en lugar de las dos que 
fueron propuestas de manera original. 

Según Killani y Ogunrombi (1984); Mutta- 
mara y Puetpaiboon (1996, 1997); Rolim (2000), 
Sperling, Chernicharo, Soares y Zerbini (2002); 
Shilton y Harrison (2003a); Shilton y Mara 
(2005); Oakley (2005); Abbas et al. (2006), y 
Cortés et al. (2013), esta condición favorece al 
flujo pistón dentro de la laguna, por tanto au- 
menta la eficiencia en la eliminación de conta- 
minantes. Los resultados del presente estudio 
concuerdan con las afirmaciones de los autores 
citados. 

En los cuadros 1 y 3 se indican las áreas 
definidas en la laguna facultativa por cada uno 


de los métodos. Las dimensiones del estanque 
facultativo con el modelo matemático son 
menores. Existe una diferencia de 594 metros 
cuadrados: 13.45%. No obstante que el área es 
menor que en el diseño tradicional, la eficiencia 
de remoción de coliformes fecales resultó ma- 
yor con el modelo de optimización. Como ya 
se mencionó, este resultado se debe a la incor- 
poración de cuatro mamparas en lugar de las 
dos propuestas de manera original. Olukanni 
y Ducoste (2011) concluyen que en la optimi- 
zación en el diseño de sistemas lagunares es 
posible reducir el costo, condicionado a que 
se incluyan las restricciones adecuadas, consi- 
derando siempre las normas de la calidad del 
agua tratada en el efluente de estos sistemas. 
Como se observa, el presente estudio coincide 
con las afirmaciones de estos autores. 

Acerca de la DBO,, los dos casos estudiados 
cumplen con la concentración de contami- 
nantes en el efluente del estanque facultativo 
que menciona la norma mexicana NOM-001- 
ECOL-96 (DOF, 1996); 43 y 48 mg/l para el 
método tradicional y la aplicación del modelo 
matemático, respectivamente. De acuerdo con 
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Cuadro 2. Informe de resultados de optimización y comparativo del diseño del sistema lagunar. 


Microsoft Excel 12.0 Answer Report 

Worksheet: [Artículo 1 solver.x1s] DISEÑO facultativa Maduración 

Report Created: 22/01/2013 16:29:42 

Result: Solver found a solution. All constraints and optimality conditions are satisfied 
Engine: Standard LSGRG Nonlinear 

Solution Time: 05 seconds 

Iterations: 6 

Subproblems: 17 

Incumbent Solutions: 3 


Objective Cell (Min) 
Cell Name Original Value Final Value 
$J$20 Costo total 5 916 989.864 5 204 275.071 


Decision Variable Cells 


Cell Name Original Value Final Value Type 
$G56 O, 25.77 22.12 Normal 
$H5$6 EE 2.00 4.00 Normal 
$156 Oi 1.88 1.00 Normal 
$J$6 Ao 4.00 2.00 Normal 
Constraints 
Cell Name Cell Value Formula Status Slack 
$H$20 Costo laguna facultativa $4 738 058.85 $H$20 >= 0 Not Binding 4 738 058.846 
$1$20 Costo laguna de maduración $466 216.23 $1$20 >= 0 Not Binding 466 216.2251 
$J$16 DBOe 41 $J$16 <= 75 Not Binding 34.08882517 
$J$16 DBOe 41 $J$16 >=0 Not Binding 40.91117483 
$1516 NF/No 1 000.00 $116 =1000 | Binding 0 
$C$16 d 0.0290 $C$16>=0 Not Binding 0.029014794 
$H56 mm. 4.00 $H$6>=1 — | Not Binding 3 
$156 O, 1.00 $156 <= 10 Not Binding 9 
$J$6 ld 2.00 $J$6>=1 | Binding 0 
$J$6 e 2.00 $J$6<=4 | Not Binding 2 
2 $156 Du 1.00 $156>=1 — | Binding 0 
2 $GS$6 O). Pao $G$6 >= 1 Not Binding 21.12037699 
E $HS6 e 4.00 $HS6<=4 — | Binding 0 
$ $GS$6 O, 22.12 $G$6 <= 30 Not Binding 7.879623014 
E 
E los resultados de los cuadros 1 y 3, se obtuvo Laguna de maduración o pulimento 
7 un costo menor con el modelo matemático de 
y $500 269.92, que representa 9.55%. Lo anterior El cuadro 3 muestra los resultados del diseño 
- es una cifra importante, considerando la dispo- de las lagunas de maduración. Con la aplica- 
Z nibilidad de recursos económicos. ción del modelo matemático propuesto resultó 


E 
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Cuadro 3. Resultados de la optimización del diseño del sistema lagunar con el sistema Risk Solver Platform. 


Datos Laguna facultativa Laguna de maduración 
Qi Ni DBOIi T O, NoE On ona 
231 10 000 000 220 11.8 22412 4 1.00 2 
Resultados de la laguna facultativa 
Xx d kb a Wan Es Qe Ne DBOe Área 
525 0.0187 0.4648 1.3301 36.70 104.09 211.90 1571 48 3819.91 
Resultados de la laguna de maduración 
XxX d kb a Wes E Qe Ne DBOe Área 
33.76 0.0290 0.4648 1.0266 35.70 8.29 210.42 1 000.00 41 295.87 
Área total Costo, laguna facultativa Costo, laguna de maduración Costo total 
4 115.78 $4 738 058.85 $466 216.23 $5 204 275.07 


menor tiempo de retención y, por tanto, un 
área menor de 39.375%. El número de mam- 
paras calculado con el método tradicional fue 
de cuatro; mientras que en el segundo diseño 
cambió a dos. No obstante lo anterior, los dos 
casos analizados cumplen perfectamente con 
los requerimientos de calidad del agua residual 
tratada que indica la Norma Oficial Mexicana, 
así como la reducción de costos, que fue de 
$212 444.86, lo anterior significa 31.30%. Como 
puede observarse, el ahorro de recursos es 
importante. La principal desventaja de estos 
sistemas de tratamiento es la gran área de te- 
rreno que ocupan (CNA e IMTA, 2007a; 2007b). 


Del sistema lagunar en su conjunto 


La reducción del tiempo de retención fue de 
4.53 días y un área menor de 16.09%. La dife- 
rencia de costos es de $712 714.79, que significa 
12.04%. 

Es prudente indicar que no es recomen- 
dable aumentar en forma indiscriminada el 
número de deflectores; según Oke y Otun 
(2001), Shilton y Harrison (2003b), Bracho et al. 
(2006) y Winfrey et al. (2010), deben tomarse en 
cuenta los costos, por lo que los autores citados 
recomiendan se realice un análisis económico 


costo-efectividad. De acuerdo con Chávez, 
Hernández y Flores (1989), se obtienen resul- 
tados óptimos que proporcionan elementos de 
juicio para la toma de decisiones. 

El cuadro 2 muestra el comparativo de 
resultados, así como las variables de decisión 
a que fue sometido el estudio. En el mismo 
cuadro, en el apartado de restricciones, se Ob- 
serva el parámetro de los coliformes fecales en 
el efluente donde se cumplió exactamente con 
la igualdad. También se puede apreciar que el 
sistema toma el número máximo permitido por 
la restricción de mamparas en la laguna faculta- 
tiva, a diferencia del tiempo de retención en la 
laguna de maduración, en el cual se selecciona 
el límite inferior de la restricción. Lo anterior se 
interpreta que cuando el número de mamparas 
aumenta en la laguna facultativa, el tiempo de 
retención de la laguna de maduración tiende a 
cero. En el cuadro 3 se observan los cambios 
realizados por el sistema de análisis Solver en 
las 20 celdas dependientes. 


Reporte de sensibilidad 
En este análisis (figura 4) se graficó la función 


objetivo y se consideró constante el número de 
mamparas determinado por el modelo mate- 
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O sis) 


$J$20 = costo total 


Figura 4. Reporte de sensibilidad de costos con variaciones de manera independiente, 


considerando sólo el tiempo de retención. 


mático de optimización en ambas lagunas: 4 
y 2, respectivamente. Al fijar como constantes 
el número de mamparas se tienen sólo tres 
variables (O,, O,, y costo total); de esta manera 
se puede representar gráficamente. Se confirma 
que la función es no lineal, ya que la relación 
entre las variables no es proporcional. También 
se infiere que el modelo es sensible al cambio 
de las variables, pues entre mayor es el tiempo 
de retención en ambas lagunas, el costo de 
construcción aumenta; así, es congruente con 
la realidad, puesto que si estos parámetros 
son mayores, aumenta el área necesaria para 
contener las aguas residuales y, por tanto, in- 
crementa el costo. 


Análisis de sensibilidad del sistema 
lagunar 


Para verificar la sensibilidad se consideró una 
variación de más menos 10% en los principales 
parámetros (Muramatsu, 2011) y se realizó una 
gráfica de tornado para observar los paráme- 


tros más sensibles al cambio. En la figura 5 se 
analiza el costo total del proyecto. 

Como se observa, la barra más ancha es el 
volumen de agua en el influente del sistema 
de tratamiento. Éste es el parámetro más sen- 
sible al cambio; es decir, a mayor volumen más 
requerimiento de área de terreno y, por ende, 
mayor costo. En orden de importancia sigue 
el tiempo de retención hidráulico, que tiene 
influencia con las dimensiones del sistema 
lagunar. 

En la figura 6 se indica que el parámetro más 
sensible que influye para la eliminación del or- 
ganismo indicador es el tiempo de retención, o 
sea, a menor tiempo de retención hidráulico es 
mayor la concentración de coliformes fecales. 
Otro parámetro importante es la temperatura: 
a menor temperatura, mayor concentración del 
organismo indicador. 

Acerca de la demanda bioquímica de oxíge- 
no en el efluente del sistema lagunar, la figura 7 
muestra que la DBO, en el influente es la barra 
más sensible, seguida por la temperatura, y 
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$J$20 
$ 4900 000.00 $5100.000.00 $ 5 300 000.00 $ 5 500 000.00 Análisis de sensibilidad, gráfica tornado 


5 4800 000.00 $5000000.00 $5 200 000.00 $ 5 400 000.00 $ 5 600 000.00 To 
| 
Parámetros 


my ap > 
ME IC 254,1 


19.908 24.332 


E laa | 
E Trial version: Limited Evaluation License dl 
g Copyright O 2011 Frontline Systems, Inc. 
E 
Figura 5. Análisis de sensibilidad de costo total. 
$1$16 
0.00 400.00 800.00 1200.00 1600.00 2000.00 En Pr zx 
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Figura 6. Análisis de sensibilidad de número de coliformes fecales. 
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Figura 7. Análisis de sensibilidad de DBO e. 
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significa que a menor temperatura, mayor con- 
centración de la materia orgánica. La eficiencia 
de remoción de la materia orgánica tiene como 
límites de temperatura superior e inferior de 
37 y 4 grados; es decir, fuera de este rango, la 
actividad de las algas que producen oxígeno 
decrece de manera importante (Rolim, 2000; 
Oakley, 2005). En orden de relevancia sigue el 
tiempo de retención hidráulico: la barra indica 
que a menor tiempo de retención mayor es la 
concentración de la materia orgánica. 

Como puede observarse, el análisis de 
sensibilidad responde de modo adecuado a la 
variación de los parámetros involucrados en 
cada concepto. 

En la figura 8 se muestran las dimensiones 
optimizadas del sistema lagunar, determinadas 
por el modelo matemático. 


Conclusiones 


De acuerdo con los objetivos planteados, se 
obtuvieron mejores resultados con la aplica- 
ción del modelo matemático de programación 
no lineal: menor costo, tiempo de retención 
hidráulico y área. 


Es importante mencionar que el presente 
análisis matemático puede aplicarse a diferen- 
tes condiciones de diseño; es decir, a cualquier 
región; pero es necesario llevar a cabo algunos 
cambios, como costos del terreno, terraplenes, 
deflectores, evaporación y temperatura, entre 
otros. 

Es recomendable, para los países en desa- 
rrollo, la aplicación de la programación no 
lineal como herramienta auxiliar alternativa en 
la optimización en el diseño de lagunas de es- 
tabilización. Se obtienen ahorros importantes y 
se cumple a la perfección con la normatividad 
en materia de descarga de aguas residuales a 
los cuerpos receptores. 

Es prudente llevar a cabo estudios adiciona- 
les, considerando la incorporación de una la- 
guna anaerobia, luego comparar los resultados 
con el fin de precisar la menor área y costo para 
la ejecución del proyecto. 
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